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Abstract. We report on a multi-band survey for very low-mass stars and brown dwarfs in the Lupus 3 cloud with 
the Wide Field Imager (WFI) at the ESO/MPG 2.2 m telescope on La Silla Observatory (Chile). Our multiband 
optical photometry is combined with available 2MASS JHK photometry to identify 19 new young stars and 3 
brown dwarf candidates as probable members of this star forming region. Our objects are mostly clustered around 
the cloud core. Stars and brown dwarfs have similar levels of Ha emission, probably a signature of accretion. One 
object, a brown dwarf candidate, exhibits a near-infrared excess, which may indicate the presence of a disk, but 
its Ha emission cannot be confirmed due to its faintness in the optical passbands. We also find two visual pairs 
of probable Lupus 3 members that may be wide binaries. 

Key words, stars: low-mass, brown dwarfs - stars: pre-main sequence - stars: formation - stars: luminosity function, 
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1. Introduction 

In the last years our understanding of the process of brown 
dwarf formation has considerably improved. Observations 
seem not to support the possibility that these objects 
form in circumstellar disks, like planets. They instead sug- 
gest that circum(sub)stellar disks, accretion and mass loss 
may be common phenomena in the substellar mass regime 
(e.g. Muench et al. ESOOTl Natta et al. I2D021 Fernandez & 
Comeron IMm Comeron et al. 120031 Mohanty et al. 120051 
Luhman et al. 120031 Scholz & Eisloffel 12003) 1 . Hence, since 
young brown dwarfs appear to have similar properties to 
low-mass pre-main sequence stars, they may be produced 
by a similar mechanism. 

Still, some questions remain open. Theoretical models 
considering the simple collapse of molecular clouds usu- 
ally fail to reproduce the observed numbers of substellar 
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objects. It is not clear how the accretion process stops be- 
fore the central object becomes massive enough to start 
the hydrogen burning. Moreover, the observed star-brown 
dwarf and brown dwarf-brown dwarf binary frequencies 
are difficult to explain with a simple extrapolation of the 
low-mass star formation models (e.g. Close et al. [2003 
Bony et al. 12005)1 . 

It has been suggested that the characteristics of the 
environment may play an important role in the forma- 
tion of brown dwarfs. For instance, Lucas & Roche (2000) 
argue that winds from the massive stars in Orion might 
have evaporated the gaseous envelopes of the nearby star 
forming cores, preventing them from growing over the hy- 
drogen burning mass limit. This would also explain why 
few brown dwarfs are known in the Taurus star forming 
region, which does not contain such massive stars (Briceno 
et al. UPQgl Martin et al. I2001|) . In Chamaeleon I, though, 
we do not find any clear relationship between the brown 
dwarfs and the intermediate mass stars in the cloud, be- 
cause they are also found near a cloud core devoid of op- 
tical objects (Lopez Marti et al. 2004, hereafter |Paper f] ). 
The possibility of formation of free-floating substellar ob- 
jects by photo-erosion of protostellar cores has been stud- 
ied by Whitworth & Zinncckcr ( 200JJ), who conclude that 
this mechanism should be very efficient in producing large 
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Fig. 1. Location of our WFI fields (Lup3-5 to 9) in the Lupus 3 cloud. The intermediate-mass stars HR 5999 and HR 6000 in 
the centre of the cloud core are seen towards the south in our field Lup3-6. The circles and triangles indicate the positions of 
the objects with and without Ha emission, respectively. The former ones are regarded as probable cloud members. 



numbers of brown dwarfs, operating over a wide range of 
conditions. However, it requires a relatively massive initial 
core to form a single low-mass object. 

Reipurth & Clarke ( 2001 ) proposed that brown dwarfs 
are stellar embryos ejected from multiple systems by dy- 
namical interactions. This hypothesis is supported by sev- 
eral numerical simulations (e.g. Bate et al.[2J)02; Delgado 
Donate et al. 2003). In a recent paper, Kroupa & Bouvier 
( 2005J) use available observational results and theoretical 
considerations to favour the ejection scenario as the prob- 
able origin of the different brown dwarf populations of 
Taurus and Orion. However, their analysis yields a num- 
ber of brown dwarfs per star that is lower than inferred by 
independent measurements of the substellar mass function 
in the Galactic field. 

To distinguish between these suggested brown dwarf 
formation scenarios it is crucial to obtain unbiased esti- 
mates of the substellar mass function from observations 
of statistically complete and homogeneous populations of 
very low- mass stars and brown dwarfs in various regions. 
Star forming regions are of particular importance since 
their youth ensures that their members have not yet suf- 
fered from strong dynamical evolution of the entire region. 



The Lupus dark cloud complex consists of five molec- 
ular clouds, named Lupus 1 to 5 (Tachihara et al. [1996), 
seen as streams or dust lanes in optical and near-infrared 
images. At least 65 T Tauri stars are known within the 
complex (e.g Schwartz ll977l Krautter ll991l Hughes et al. 
1994|, most of them concentrated in the Lupus 3 cloud. 
Two A-type stars, HR 5999 and HR 6000, are seen to- 
wards the centre of this cloud. Their membership of the 
region is controversial, though, because the measured dis- 
tance to these stars is in disagreement with other distance 
estimates to the complex (e.g. Crawford 2000). 

A peculiarity of Lupus is that the distribution of spec- 
tral types for the stars in these clouds is dominated by 
midM-type objects, in contrast to other star forming re- 
gions like Taurus- Auriga, which contain many late K-type 
T Tauri stars, as noted by many authors (e.g. Appenzeller 
et al. IT5%3l Hughes et al. IT5M1 Wichmann et al. IlW|) . 
For this reason, Lupus appears as a very promising region 
to search for young brown dwarfs. Recently, Comeron et 
al. H2003|l reported a bona-fide brown dwarf, Par Lup3- 
1, from a small Ha survey around the stars HR 5999 
and HR 6000. A recent survey in the near-infrared by 
Nakajima et al. i2000[) has also provided some brown 
dwarf candidates. 
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Table 1. Log of WFI observations in Lupus 3 



Date Field a (2000) 8 (2000) 



28 May 1999 


Lup3-6 


16 h 


08 m 


36.5 s 


-38° 57' 


36'.' 7 


30 May 1999 


Lup3-7 


16 h 


08 m 


26.3 s 


-39° 28' 


47'.' 1 


30 May 1999 


Lup3-8 


16 h 


ll m 


01.2 s 


-38° 56' 


27'.' 9 


3 Jun 1999 


Lup3-5 


16 h 


05 m 


40.6 s 


-39° 22' 


32" 2 


3 Jun 1999 


Lup3-9 


16 h 




07.8 s 


-39° 27' 


37'.' 6 



Table 2. Exposure times (in seconds) for the images taken in 
each selected filter. 



Filter 


Short exp. 


Interm. exp. 


Long exp. 


Rc/162 


5 


60 


600 


Ic/lwp 


5 


30 


600 


Halpha/7 


15 




600 


855/20 


16 


300 


600 


915/28 


8 


100 


600 



In this paper we present the results from our multi- 
band survey for brown dwarfs in the Lupus 3 cloud. The 
method followed is the same as in our previous survey in 
Chamaeleon I QPaper l| . The current paper is structured 
as follows: Our observations are described in Sect.El The 
candidate selection and brown dwarf identification are ex- 
plained in Sect.|3 In Sect. 0] we analyse the observed spa- 
tial distribution of our objects and the frequency of bina- 
ries. Their Ha emission properties are discussed in Sect. [5] 
Finally, Sect. El summarises our conclusions. 

2. Observations 

An area of about 1.6 D was observed in the Lupus 3 cloud 
with the Wide Field Imager (WFI) at the ESO/MPG 
2.2 m telescope on La Silla Observatory (Chile). The sur- 
veyed area corresponds to five WFI fields, each of 34' x 33' 
in size (see Fig. QJ. Four of the fields (Lup3-6, 7, 8 and 
9) define a rectangle centred on the cloud core, seen as 
an elongated dust lane in our images. The two A-stars 
HR 5999 and HR 6000 are placed in the southern part of 
the field Lup3-6. The Lup3-5 field lies towards the south- 
west of this area, also following the orientation of the core. 
Table H contains the coordinates of the WFI field centres 
and the date of observation for each field. The observations 
were carried out in the same three nights as our analogous 
survey in Chamaeleon I ( |Paper l\ . 

We observed in two broad-band filters, R and I, in a 
narrow-band filter centred in the Ha emission line (Ha/7), 
and in two medium-band filters, M855 (855/20) and M915 
(915/28). The last two allow for a photometric spectral 
type classification (see Sect. below). For each field and 
filter, except for Ha, three images with different exposure 
times were taken (see TableEJ- In the case of Ha, we took 
only two different exposures per field. 



Data reduction, object search, and photometry were 
performed within the IRAF environment 1 in the way ex- 
plained in |Paper 1| After the standard reduction (bias 
subtraction and flatfield division) , all images were divided 
through an illumination mask to correct for an irregular 
illumination pattern. For the reddest filters (I, M855, and 
M915), it was also necessary to subtract a fringe mask 
to correct for the strong fringing pattern. For each filter, 
both masks were created by combination of the science 
exposures as explained by Scholz & Eisloffel (20T)4^. 2 

We ran SExtractor (Bertin & Arnouts 1996) to pro- 
duce an object catalogue of each surveyed field. For the 
photometry, the DAOPHOT package was used (Stetson 
1987). PSF photometry was performed in order to also 
measure close pairs of objects. The broad-band photom- 
etry was calibrated with Landolt l|1992f) standard stars. 
The extinction coefficients from the fit of these stars for 
the R and I bands were also used to perform a correc- 
tion for atmospheric extinction of the Ha and the M855 
and M915 photometry, respectively. No set of standard 
stars is available for our narrow and medium-band filters. 
For a better understanding of the emission properties, the 
Ha instrumental magnitudes were then shifted so that the 
main locus of the objects in each surveyed field corre- 
sponds to Ha-R= in a (Ha, Ha-R) colour-magnitude 
diagram (see Sect.[3J|. The mean error in this shift is of less 
than 0.01 mag for the fields Lup3-6, Lup3-7 and Lup3-8, 
and of about 0.05 mag in the fields Lup3-5 and Lup3-9. 
This value is comparable to the mean photometric error 
(see below) . We estimate that our survey is complete down 
to R~20 mag and I~19 mag. 

As discussed by Alcala et al. ( 2002J , the WFI chip-to- 
chip photometric variations can be as high as 3% in the 
broad-band filters and 5% in medium-band filters. This is 
the main error source in our photometric calibration, since 
only a global fit for all the chips could be performed with 
the Landolt standard stars. Moreover, a systematic offset 
was found between images with different exposure times, 
in the sense of the resulting luminosities being brighter 
(typically by about 0.05 mag) for longer exposure times. 
The origin of this offset is so far unknown. To minimise 
this systematic error, we determined the offsets between 
the three different integration times, and then shifted the 
short and long integration times to the system of the 
medium integration. 

The errors in the completeness range are in general not 
larger than 5%. A more detailed discussion of the calibra- 
tion procedure and the sources of the photometric errors 
can be found in |Paper 1| 



1 IRAF is distributed by the National Optical Astronomy 
Observatory (NOAO), operated by the Association of 
Universities for Research in Astronomy, Inc. under contract 
to the US National Science Foundation. 

2 An alternative method to correct the WFI illumination gra- 
dient using a second-order polynomial has been proposed by 
Koch et al. JSPOIll . 
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Due to problems with the shutter, several Lup3-7 and 
Lup3-9 images (including the long and intermediate I- 
band exposures, necessary for the object selection) exhibit 
some darker regions at the field edges. Objects in these re- 
gions will not be considered in the following analysis. 

3. Brown dwarf candidates 

3.1. Candidate selection 

In Paper 1, candidate members of Chamaeleon I were 
selected around an empirical isochrone defined by the 
previously known brown dwarfs and very low-mass stars 
in the cloud in a (I, R-I) colour-magnitude diagram. In 
Lupus 3, however, there is no sample of known brown 
dwarfs that could be used for this purpose. Hughes et al. 
(1994) provide broad-band RIJHK photometry of a sam- 
ple of known T Tauri stars in the cloud, but most of these 
bright objects are saturated in our images. Therefore, 
we had to rely on the theoretical models in this case. 
Candidate members of Lupus 3 were thus selected around 
an isochrone from the models of Chabrier et al. ( 2000 ) . To 
do this, we needed an estimation of both the extinction 
and the distance to this star forming region. 

Hughes et al. l|1994f) estimated the individual extinc- 
tion values towards the Lupus stars in their study. We used 
these data to derive an average extinction of Av=0.86iQ g^ 
in our surveyed region. With this value, and using the fol- 
lowing relations: 

Ai/A v = 0.482 (1) 
Ai = 1.812-[(R-I)-(R-I) ], (2) 

derived from the extinction law of Rieke & Lebofsky 
(1985 ), an average colour excess of 

(R - I) -(R- 1)0 = 0.231°;^ (3) 

was computed for the Lupus objects. It must be noted 
here that Rieke & Lebofsky's work uses the Johnson R 
and I filters, while the system used in this paper is closer 
to the Cousins one. However, tests using the extinction 
law of Mathis (1990) showed that the sample of selected 
objects would not be affected by the use of one or the 
other extinction law. 

There has been some debate about the actual distance 
to the Lupus clouds. Values between 130 and 220 pc can be 
found in the literature (see e.g. Crawford 2000]for a discus- 
sion). Since our extinction estimate is based on the work 
by Hughes et al. (199||), we will consider, for consistency, 
the distance given by these authors, 140±20 pc. Within 
the quoted errors, this value is consistent with most of 
the more recent measurements (see e.g. Crawford 2000 
Comeron et al. 2003). However, as discussed later in this 
paper, this might be an underestimation of the true dis- 
tance to Lupus 3 (see Sect. l3.2l belowL 

Fig. |21 shows the (I, R-I) colour-magnitude diagram 
for our surveyed field in Lupus 3. The solid lines are a 




R-[ 



Fig. 2. (I, R-I) colour-magnitude diagram for our surveyed 
field in Lupus 3. For clarity, only a random selection of the 
background objects has been plotted. Brown dwarf candidates 
are selected around two isochrones from the models of Chabrier 
et al. 1)2000 ) corresponding to the ages of 1 and 5 Myr (thin 
and thick solid line, respectively). Objects with and without de- 
tected Ha emission are marked with solid and open circles, re- 
spectively. The dashed lines indicate the range in colour around 
the isochrones taken as limit for the candidate selection, while 
the dotted line is the empirical isochrone used for candidate 
selection in Chamaeleon I. The crosses indicate the average er- 
rors. Note that objects fainter than 1=20.0 mag (correspond- 
ing to planetary masses) are not considered in the candidate 
selection. The theoretical models place the stellar/substellar 
boundary at I~ 14 mag for the distance, extinction, and age 
interval considered. 

1 Myr and a 5 Myr isochrone from the Chabrier et al. 
(2000) models, shifted to the estimated distance and ex- 
tinction of the cloud. These isochrones span the age range 
of the Lupus stars. The dashed lines indicate the position 
of the isochrones according to the upper and lower limits 
for the extinction derived in the above calculation. A large 
number of objects are found between these two lines. Note 
that there are very few objects to the right of this area; 
the density of objects is also clearly diminished left of our 
search area with respect to the main locus of objects seen 
in the field. Hence, our extinction estimate appears as a 
reasonable approximation. 

In the Lup3-5 field, objects on top of the third CCD 
chip have bad flux measurements in the long R expo- 
sure due to a badly subtracted background. The conse- 
quence was that these objects had too faint R magni- 
tudes and thus appeared very red in the (I, R-I) diagram 
and very blue in the (Ha, Ha-R) diagram. Visual inspec- 
tion showed, however, that these objects were similarly 
bright in R and I, and thus very unlikely cloud mem- 
bers. Therefore, they were rejected from our candidate 
list, and are not included in the colour-magnitude dia- 
grams of Figs. [21 El and 0] Excluding these objects, we 
have 85 candidates. 

We can compare the diagram in Fig. |21 with the (I, R- 
I) colour-magnitude diagram for Chamaeleon I (Fig. 2 in 
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Fig. 3. (Ha, Ha-R) colour-magnitude diagram for our sur- 
veyed field in Lupus 3. For clarity, only a random selection of 
the background objects has been plotted. Our candidates with 
and without Ha emission are plotted as solid and open circles, 
respectively. Ha emitters are placed in the left part of the dia- 
gram, having negative Ha-R colour index. (Note that objects 
not visible in our Ha images could not be plotted.) 

|Paper 1\ . The dotted line in Fig. [5] indicates the position 
of the empirical isochrone used for candidate selection in 
|Paper l| Not unexpectedly, given that both regions are 
at similar distances, this line is almost coincident with 
the 1 Myr isochrone, except for the faintest objects, for 
which the R-I colour saturates according to the theoretical 
models. However, as discussed in detail in |Paper l| the 
uncertainties in the colours for the faintest objects are 
large, due both to the photometric errors and to the lack of 
standard stars that are red enough for a good photometry 
calibration in this range. This fact may have implications 
for the contamination with background objects among our 
faintest candidates, as discussed below. 

All the found objects are fainter than I~14 mag. 
According to the Chabrier et al. (2000) models, this 
magnitude roughly corresponds to the stellar/substellar 
boundary (M~0.075 M Q ) at the assumed age (~2 Myr, 
see e.g. Wichmann et al. I1997fl and distance of Lupus 3 
(140 pc). The effective temperature given by these mod- 
els for objects with this mass at this age is around 2900- 
3000 K, corresponding to spectral types around M6 ac- 
cording to the scale of Luhman ( 1999). Thus, if these can- 
didates are indeed members of the star forming region, 
they must lie close or beyond the star/brown dwarf tran- 
sition. 

Since Lupus 3 is a nearby dark cloud, we do not ex- 
pect, in principle, a high contamination from foreground 
and background objects. Simulations of the Galaxy pop- 
ulation using the Besangon models (Robin & Creze 1986 
Robin et al. l2003f) towards the direction of the Lupus com- 
plex yield between 6 and 15 objects in the region of in- 
terest of the (I, R-I) colour-magnitude diagram. Taking 
the most conservative result, we would expect that only 
around 15% of our candidates should be non-members of 



the star forming region. However, the colour uncertainties 
for our faintest objects mentioned above very probably in- 
crease this percentage at the lower magnitude limit of our 
survey. Given the large errors for the faintest objects, we 
stopped our candidate selection at I~ 20 mag in Fig. 
Since the Chabrier et al. (2000) models place the plan- 
etary mass limit at around this value for our estimated 
age, distance, and extinction, we do not expect to have 
planetary mass objects among our candidates. 

To confirm the youth of our objects, we tested their Hcv 
emission. As shown in |Paper"Tl there is a relation between 
the Ha-R colour index and the logarithm of the equivalent 
width of the Hcv emission line, logEW(Ha): In general, the 
bluer (i.e., more negative) the Ha-R colour, the stronger 
the Ha emission. Observations of previously known brown 
dwarf candidates in Chamaeleon I showed that objects 
with Ha equivalent widths of llA already present Ha-R 
colours bluer than the main locus of objects in the field 
(see |Paper l| for details). 

Fig. El shows the (Ha, Ha-R) colour-magnitude dia- 
gram for our Lupus survey. Many of the brightest objects 
exhibit clear Ha emission, i.e. negative values of the Ha-R 
colour index. Some objects are clear non-emitters, having 
positive values of the Ha-R colour index. A number of 
the faintest candidates also seem to have Ha emission, 
but this is very uncertain due to the large photometric 
errors; therefore, except for an object with an apparent 
near-infrared excess (see Sect. 13.31 belowL they were not 
further considered. Moreover, twenty-two of the faintest 
objects are not detected in our Ha images; hence, they 
do not appear in Fig. [3] Since nothing can be said about 
their emitting properties, these 22 faintest objects were 
also excluded from our study. 

For the brightest objects in Fig.OI we applied the same 
criterion as Lamm et al. 1|2004|) . retaining the objects with 

A (Ha - R) = (Ha - R) bjcct - (Ha - R) locus < -0.1 (4) 

as probable cloud members. We discarded all the faintest 
candidates (Ha >20.5 mag), because the large errors at 
the lower part of the colour-magnitude diagram do not al- 
low us to discriminate between emitters and non-emitters. 

Our final sample contains 22 stars and brown dwarf 
candidates. Note that more than 65% of the objects de- 
tected in our Ha images do not seem to have Ha emission 
according to this criterion. Such a low fraction of objects 
with Ha emission is not surprising. In NGC 2264, which 
has an age (1-3 Myr) very similar to Lupus 3, Lamm et 
al. (2004) found only 18% of strong Ha emitters among 
their low-mass stars with masses M < 0.25M. On the 
other hand, some objects might have Ha emission levels 
too low to be reliably detected with our data. Moreover, 
not all young objects show Ha emission. The scarcity of 
Ha emitters among very low-mass stars in other star form- 
ing regions suggests that many of the remaining objects 
from our initial candidate sample belong to the star form- 
ing region. Hence, future observations are likely to increase 
the number of very low-mass members in Lupus 3. 
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M855-M915 

Fig. 4. (M915, M855-M915) colour-magnitude diagram for 
our surveyed field in the Lupus 3 cloud. Objects with and 
without Ha emission are marked with solid and open circles, 
respectively. Our scale for the identification of the spectral type 
is also indicated. The faintest objects in the diagram are re- 
garded as possible L-type objects. 

3.2. Spectral types 

Our medium-band filter M855 lies in a spectral region in- 
cluding the TiO and VO absorption bands used to clas- 
sify M-type objects (see Kirkpatrick et al. I1991|) . The fil- 
ter M915 covers a pseudocontinuum region in these late- 
type objects. In |Paper"T| we showed that there is a rela- 
tion between the M855-M915 colour index and the M- 
subspectral type using synthetic photometry from the 
models of Baraffe et al. (|1998» . In the range M4-M9, a 
linear fit is a good approximation of this relation. The cor- 
responding correlation coefficients were determined using 
photometry of known objects of late spectral type. On the 
other hand, for spectral types earlier than M4 the M855- 
M915 colour saturates, because the TiO and VO features 
are not prominent or not present at all in the spectra. A 
second correlation, of opposite sign, was found for spec- 
tral types L0-L2, due to the progressive disappearance of 
the TiO and VO from the spectra. In this case, however, 
the uncertainties are larger than for late M-type objects, 
given that extincted early M-type objects would be placed 
in the same region of the colour-magnitude diagram. 

This method enabled us to derive spectral types for our 
brown dwarf candidates from their position in a (M9I5, 
M855-M915) colour-magnitude diagram. We estimated an 
error of 1-2 subclasses in this calibration. A detailed dis- 
cussion of this photometric spectral type classification, its 
error sources, and its limitations can be found in |Paper"l"1 

Fig. H shows the corresponding (M9I5, M855-M915) 
colour-magnitude diagram. The vertical lines indicate our 
scale for the spectral type determination. It is evident 
from this figure that most of our objects have spectral 
types M4-M5, thus lying above, but close to, the substel- 
lar boundary. A few objects have spectral types M6 or 
later, being brown dwarfs or brown dwarf candidates. A 



"turning back" towards the left, which could be indicative 
of the beginning of the L spectral sequence, is seen for 
the faintest objects. This change of tendency happens at 
a similar M91 5 magnitude as in the case of Chamaeleon I 
(Fig. 8 in Paper I), as expected if both clouds have similar 
age, distance, and average extinction. 

About one half of the fainter candidates are found at 
the lower left part of the diagram, in the main locus of 
the objects seen in the surveyed field. In this group, only 
Lup 654 and Lup 706 seem to have Ha emission. The other 
two or three objects in the lower part of the diagram do 
not show clear Ha emission in Fig. [31 Except for Lup 642, 
which shows a near-infrared excess and apparently no ex- 
tinction (see below) , they could not be confirmed as cloud 
members, and are thus rejected. 

The WFI photometry results, spectral types, and clas- 
sification for our 22 new Lupus 3 members are summarized 
in Table |3| As in |Paper"l"l all objects with spectral types 
earlier than M5.5 (16 objects) are classified as very low- 
mass stars. Four objects lie in the transition from stars to 
brown dwarfs with spectral types between M6 and M6.5. 
Three objects Lup 642, Lup 654, and Lup 706, have esti- 
mated spectral types later than M6.5. Although they are 
considered brown dwarf candidates in the subsequent dis- 
cussion, it must be kept in mind that their spectral types 
(all later than M9) are very uncertain. 

Note that the largest group of objects corresponds to 
spectral type M5, as expected from the characteristics of 
the Lupus region. As mentioned in Section^ the Lupus 3 
cloud is particularly rich in mid M-type objects. It is re- 
markable, though, that the new members in Tableware in 
general fainter than expected for their estimated spectral 
types: As stated in the previous section, with the assumed 
distance and extinction, an object with I~14 mag should 
be of spectral type M6 according to the theoretical mod- 
els. However, for such objects (see e.g. Lup 915 in Tablc|2J 
we derive a spectral type M4, about 2 subclasses earlier 
than the theoretical prediction, although still marginally 
coincident with the expectations (we estimate an error of 
about 2 subclasses in our classification). 

The discrepancy between the luminosities predicted by 
the theoretical models and our results cannot be explained 
by the photometric errors alone. It may be a consequence 
of the underestimation of the average extinction in the 
cloud (see Sect.EUl. To obt ain the measured magnitudes, 
Ay should be about 4.5 mag, a value not unusual in star 
forming regions. This would also affect our derived spec- 
tral types, because the objects would appear redder and 
fainter, and hence of later spectral type, in our colour- 
magnitude diagram (see discussion in |Paper"H ). However, 
most of our objects do not seem to be highly extincted (see 
Sect. 13. 3)l . It may also be that the distance to Lupus 3 has 
been underestimated. At 200 pc and with the extinction 
estimate from Sect. [31 a M4 object would have I~13 mag, 
which is closer to our measured values. Indeed, compari- 
son of the predictions of the Lyon models (Baraffe et al. 
1998) with the photometry and spectral types of the mid 
M-type stars studied by Hughes et al. 1|1994J> . whose dis- 



Lopez Marti' et al.: Very low-mass members of the Lupus 3 cloud 



7 



tance value has been taken in this work, seems to indicate 
that, if the extinctions given by these authors are correct, 
the objects are too faint for their spectral types. For in- 
stance, Sz 112, a M4 star with an estimated extinction of 
A v =0.85 mag, very close to our assumed average value, 
has according to Hughes et al. Q1994[) an I-magnitude of 
13. This also hints to a wrong distance estimate. 



A further possibility would be that these objects ap- 
pear underluminous for their spectral types due to sur- 
rounding dust or accretion. Such objects have been re- 
ported in other star forming regions (e.g. Fernandez & 
Comeron 2001; Luhman et al. 2003). Although it seems 
unlikely that most of our Lupus members are affected by 
such underluminousity, only reported for a minority of ob- 
jects so far, the possibility cannot be completely excluded 
without a safer distance and age determination. 



There are not many previous optical surveys in 
Lupus 3 that are deep enough to be compared with our 
study. We recovered one object from the list provided by 
Hughes et al. 1|1994[> . Sz 113 (Lup 609s). Our derived spec- 
tral type for this object (M5) is later but still consistent, 
within our estimated errors, with the classification of these 
authors (M4). On the other hand, Comeron et al. (2003) 
estimate an even later spectral type (M6) for the same 
object. They suggest that the discrepancy between their 
result and that of Hughes et al. is due to the different 
spectral ranges used to perform the classification. In view 
of this, we decided to keep our derived spectral type for 
this object. The classifications by other authors are given 
in a footnote to Table |3J Other objects from Hughes et 
al. I|1994fl seen in our images are too bright and hence 
saturated. 



Comeron et al. (|2003f) reported on four new Lupus 3 
members identified in a spectroscopic survey of a small 
area around HR 5999 and HR 6000. Of them, one (Par- 
Lup3-2) is not present in our images because it falls in 
an interchip gap. The other three have R-I colours that 
place them at the edges of our selection band in Fig. [2] 
While Par-Lup3-4 has a relatively blue R-I colour that 
may be caused by its very high Ha emission, Par-Lup3-1 
and Par-Lup3-3 are placed to the red of our objects in the 
(I, R-I) colour-magnitude diagram. Since these objects are 
all placed in a small area around the stars HR 5999 and 
HR 6000, in a region of high extinction in the Lupus 3 
cloud, this might indicate that we are missing highly ex- 
tincted objects with our selection criteria. However, as 
commented in Section l3.ll there are not many objects to- 
wards the right of our selection band, so this loss should 
not be very significant. Moreover, objects with very red 
R-I colours (>3 mag) are more likely to be background 
stars seen through the dark cloud. 



Av-10 

o Av=5 

2 - ° O Av=0.86 - 




- Lup642 

O 

0.4 0.8 1.2 

H-K 

Fig. 5. Colour-colour diagram for our candidates with near- 
infrared photometry. Objects with and without Ha emission 
are marked with solid and open symbols, respectively. The solid 
curves indicate the locus of dereddened dwarfs and giants ac- 
cording to Bessell & Brett |1988|. The dotted lines indicate 
the direction of the reddening vector up to Av ~10 mag. The 
other lines indicate the position of a 1 Myr isochrone from the 
models of Baraffe et al. 0.998 ) for different extinction values. 

Table 4. 2MASS photometry for our candidate members of 
Lupus 3 



Name 


J 


H 


K 


2MASS Name 




mag 


mag 


mag 




Lup 604s 


12.145 


11.452 


11.071 


16080017-3902595 


Lup 608s 


12.202 


11.642 


11.358 


16090850-3903430 


Lup 609s 


12.465 


11.724 


11.257 


16085780-3902227 


Lup 605 


14.207 


13.569 


13.201 


16071400-3852375 


Lup 617 


12.870 


12.257 


11.965 


16084817-3909197 


Lup 642 


15.381 


15.536 


14.162 


16090155-3905055 


Lup 710s 


13.548 


13.044 


12.661 


16091713-3927096 


Lup 713s 


13.240 


12.566 


12.130 


16073773-3921388 


Lup 706 


15.174 


14.240 


13.829 


16083733-3923109 


Lup 714 


12.576 


11.922 


11.590 


16075888-3924347 


Lup 802s 


13.252 


12.561 


12.254 


16115116-3851047 


Lup 810s 


12.837 


11.683 


11.263 


16095449-3912035 


Lup 818s 


13.006 


12.392 


11.991 


16095628-3859518 


Lup 831s 


12.507 


11.192 


10.689 


16113859-3908273 



3.3. Infrared detections 

We made use of the database of the 2MASS survey 3 
to look for near-infrared counterparts of the 102 objects 
(with and without Ha emission) in our initial candidate 
list. In this way, we intended not only to further con- 
firm the youth of our 36 Ha emitters, but also to even- 
tually identify new candidate cloud members. We found 
40 positive detections, of which 27 are objects without de- 
tectable Ha emission from their Ha-R colours. All 2MASS 
sources are placed at distances not larger than about 1" 
from the given positions for our objects, hence the cross- 
identifications are very reliable in all cases. 

Table 0] summarizes the J UK photometry of the new 
Lupus 3 members from Table that were detected by 
2MASS. The photometric errors are in general around 

3 Available online at the URL of the NASA-IPAC Infrared 
Science Archive (IRSA): http://irsa.ipac.caltech.edu/ 
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Table 3. Candidate low-mass members of Lupus 3 ab 



Name 


a (2000) 


8 (2000) 


R 


I 


Ha 


M855 M915 SpT 1 " Classification 




hh mm ss.s 


ddd mm ss.s 


mag 


mag 


mag 


mag mag 



Lup 


504 


16 


06 


47 





-39 


16 


15 


8 


16.77 


14 


94 


16.38 


18 


54 


18 


00 


M4 


star 


Lup 


604s 


16 


08 


00 


2 


-39 


02 


59 


7 


16.51 


14 


36 


16.09 


18 


34 


17 


65 


M5.5 


star 


Lup 


608s 


16 


09 


08 


5 


-39 


03 


43 


8 


16.30 


14 


24 


15.99 


18 


21 


17 


57 


M5 


star 


Lup 


609s d 


16 


OS 


57 


8 


-39 


02 


23 


6 


17.05 


14 


88 


16.34 


18 


81 


18 


14 


M5° 


star 


Lup 


605 


16 


07 


14 





-38 


52 


37 


9 


18.96 


16 


46 


18.45 


20 


41 


19 


58 


M6.5 


trans, obj. 


Lup 


007 


16 


08 


28 


1 


-39 


13 


09 


6 


18.39 


16 


10 


18.06 


19 


94 


19 


36 


M5 


star 


Lup 


617 


16 


08 


48 


2 


-39 


09 


20 


1 


17.27 


15 


02 


16.94 


18 


97 


18 


22 


M6 


trans, obj. 


Lup 


642* 


16 


09 


01 


5 


-39 


05 


06 


1 


22.6: 


19 


.9; 


21.7: 


23 


.0: 


23 


0; 


L2 


BD cand. 


Lup 


648 


16 


09 


48 


6 


-39 


11 


17 


6 


17.11 


15 


28 


16.69 


19 


25 


18 


64 


M5 g 


star 


Lup 


650 


16 


09 


49 


8 


-38 


49 


04 


5 


19.84 


17 


62 


19.5; 


21 


36 


20 


85 


M4 


star 


Lup 


652 


16 


07 


09 


5 


-38 


41 


30 


3 


20.7; 


18 


36 


20.1: 


22 


•i; 


21 


54 


M4.5 


star 


Lup 


654 


16 


07 


23 


4 


-39 


05 


13 


2 


19.92 


17 


25 


19.5; 


22 


.5; 


22 


3; 


LI 


BD cand. 


Lup 


710s 


16 


09 


17 


1 


-39 


27 


09 


4 


17.98 


15 


84 


17.60 


19 


80 


19 


15 


M5 


star 


Lup 


713s 


16 


07 


37 


7 


-39 


21 


38 


8 


17.61 


15 


66 


16.28 


19 


69 


18 


90 


M0 


trans, obj. 


Lup 


706 


16 


08 


37 


3 


-39 


23 


10 


8 


20.7; 


18 


27 


19.1; 


22 


.5; 


21 


24 


LO 


BD cand. 


Lup 


707 


16 


08 


28 


1 


-39 


13 


09 


7 


18.11 


16 


09 


17.60 


20 


00 


19 


39 


M5 


star 


Lup 


714 


16 


07 


58 


9 


-39 


24 


34 


9 


16.78 


14 


81 


16.51 


18 


80 


18 


14 


M5 


star 


Lup 


802s 


16 


11 


51 


2 


-38 


51 


04 


2 


17.11 


15 


29 


16.78 


19 


30 


18 


84 


M4 


star 


Lup 


810s 


16 


09 


54 


6 


-39 


12 


03 


4 


16.60 


14 


84 


16.42 


19 


20 


18 


81 


<M4 


star 


Lup 


818s 


16 


09 


56 


3 


-38 


59 


52 


1 


17.42 


15 


25 


16.60 


19 


34 


18 


57 


M6 


trans, obj. 


Lup 


831s 


16 


11 


38 


6 


-39 


08 


27 


1 


17.25 


15 


32 


17.14 


19 


23 


18 


82 


<M4 


star 


Lup 


914 


16 


10 


16 


1 


-39 


37 


53 


4 


17.87 


15 


81 


17.64 


20 


79 


20 


62 


<M4 


star 


Lup 


915 


16 


10 


54 


1 


-39 


40 


07 





15.93 


14 


05 


15.72 


18 


70 


18 


24 


M4 


star 



Ha, M855, and M915 magnitudes arc instrumental magnitudes. 
Photometric errors: blank: 0.05 mag; semicolon: 0.1 mag; colon: 0.2 mag 
Errors in the spectral types: M4-M6: 2 subclasses; M6.5-M9: 1 subclass 
This object is identified with Sz 113 (Schwartz 1 9771 

Hughes et al. 119941 and Comeron et al. I20U3I give a spectral type of M4 and M6, respectively. 
Near-infrared selected (see text). 
Ambiguous spectral type (see text). 



0.02-0.03 in all three bands. In total, about 50% of the 
stars and all the transition objects from Table are de- 
tected by 2MASS. Two of the brown dwarf candidates are 
also found. 

Above I~19 mag, more than 70% of the objects were 
detected by 2MASS, the missing objects in this range be- 
ing all placed in the outer parts of the cloud. This is not 
unexpected, because in those regions the K magnitude for 
a given I magnitude should be fainter due to lower red- 
dening. Most of the WFI candidates not detected in the 
infrared belong to the faintest optical sources; hence, they 
probably lie below the 2MASS detection threshold. 

Fig. shows the (J — H, H — K) colour-colour dia- 
gram for our 2MASS objects. Objects with and without 
Ha emission are marked with different symbols. The solid 
lines indicate the position of the main sequence stars ac- 
cording to Bessell & Brett II1988|) . The positions of a 1 Myr 
isochrone from the Baraffe et al. (1998) models for several 
values of extinction are also indicated. 

A degeneracy in the near-infrared colours exists be- 
tween main sequence late-type dwarfs and giants and pre- 
main sequence low-mass stars. On the other hand, no mat- 
ter whether they are true Lupus members or not, the ob- 
jects in Fig. are not equally extincted. All the objects 
without Ha emission have moderate to high values of Ay- 
(>5 mag). Most of them are found along the reddening line 
for giants (upper dotted line); hence, given the lack of Ha 
emission, they are probably obscured stars seen through 
or at the edges of the dark cloud. In contrast, most of 
the Ha emitters are found at low or zero values of extinc- 



tion, consistent with our selection criterion from Fig. 
Three objects with Ha emission are placed at values of 
Ay between 2 and 5 mag. 

Only one object, Lup 642, seems to have a near- 
infrared excess according to the colour-colour diagram in 
Fig- El This detection suggests its belonging to the Lupus 3 
cloud, despite not having clear Ha emission. However, this 
object lies very near a previously known Lupus T Tauri 
star (see Sect. FQJ, also detected by 2MASS, and its 
near-infrared photometry might be contaminated by this 
brighter source. 

No other object from our initial sample shows excess 
emission in the near-infrared. Such an excess is commonly 
detected in young objects and has its origin in a cir- 
cum(sub)stellar disk. It must be noted, however, that the 
non-detection of a near-infrared excess does not automati- 
cally exclude an object from having a disk. We will discuss 
this fact later in Sect.03 

4. Spatial distribution 

4.1. Observed distribution 

Fig.fflshows the positions of our new very low-mass mem- 
bers of the Lupus 3 cloud (dots) and of the objects without 
Ha emission from our initial selection (triangles). Most 
of the objects are found within or at the boundaries of 
the cloud core (seen as an elongated dust lane in the im- 
age), particularly near the apparent location of the two 
intermediate-mass stars HR 5999 and HR 6000. 
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Fig. 6. Distribution of the low-mass Lupus members with dis- 
tance to the pair of intermediate mass stars HR 5999 and 
HR 6000. The blank histograms show the distribution of the 
very low-mass stars; the hashed histograms, of the brown 
dwarfs and faint brown dwarf candidates. The common his- 
togram for all these objects is showed in dotted lines. 

In analogy to the analysis performed for Chamaeleon I, 
we have counted the number of objects at different dis- 
tance bins from the intermediate-mass stars; the his- 
tograms derived from these counts are shown in Fig. 
Two maxima are found at 5-10' and 15-20' from the inter- 
mediate mass stars, corresponding to the cluster of very 
low-mass objects seen around HR 5999 and HR 6000 in 
Fig. n] The three brown dwarf candidates are all placed 
within this distance. Further away, the numbers of stars 
decrease. We remark that the bins include all objects at 
a given distance in all directions. Inspection of Fig. ^ 
shows that, in fact, the objects further than about 15- 
20' are not concentred in one particular area. At shorter 
distances, however, most of the objects are seen along the 
dark ridge, a region of higher density in the cloud (see 
e.g. Tachihara et al. 1996). Note also that there might 
be some incompleteness in the first bin (distances 0-5'), 
because very close faint objects might be blended by the 
brighter intermediate-mass stars. 

Although some objects are found far from the cloud 
core, it is clear from Figs. and |H| that most of the 
Lupus members, and particularly the brown dwarf can- 
didates, are still very close to their birth site. This dis- 
tribution does not seem to support the ejection model of 
Reipurth & Clarke (2001) unless, as argued by Kroupa & 
Bouvier ( 2003 ) , dynamical interactions between members 
are able to keep the ejected brown dwarfs within the clus- 
ter. Kroupa & Bouvier (2003) successfully explain why 
a large number of these objects are found in the star- 
rich Orion clusters, but not in the more sparse Taurus- 
Auriga star forming region. However, Lupus (and also 
Chamaeleon I) is an intermediate case between these two. 
Detailed calculations for this region should be performed 
before stating whether dynamical interactions could have 
retained the ejected brown dwarfs also in this cloud. 



The association with the intermediate mass stars 
might be accidental, since their belonging to the cloud 
is not clear: Several distance measurements yield a value 
of about 190-220 pc for these stars (see Wichmann et al. 
1998; Crawford 20® and references therein), inconsistent 
with other estimations of the distance to the Lupus clouds 
(between 140 and 160 pc). On the other hand, in Pa per 1| 
no clear relationship could be established between the 
very low-mass objects and the intermediate mass stars in 
Chamaeleon I; the former were found near the dark cloud 
cores, independently of the presence of higher mass stars 
in these cores. 

4.2. Wide binaries 

To check the existence of multiple systems among our ob- 
jects, we first derived the average surface density in our 
surveyed area. Then we computed the minimum distance 
at which an object should have a neighbour if the distri- 
bution were completely uniform. In this way, we estimated 
a maximum separation of about 5'.' 5 (about 770 AU for a 
distance of 140 pc) to look for binary candidates. 

Despite their clustering around the cloud cores, few 
visual pairs are found among our probable Lupus mem- 
bers. Only two objects, the transition object Lup 818s 
(M6) and the faint brown dwarf candidate Lup 642, are 
seen near a previously known Lupus star: The former lies 
at about 4" (~560 AU) of Sz 119, while the latter is 
found at about 6" (-840 AU) of Sz 114. The visual bi- 
nary frequency of Lupus 3 is thus very similar to that of 
Chamaeleon I, where also only a few possible systems were 
seen ( |Paper l\ . 

Only precise proper motion and radial velocity mea- 
surements showing common space motion may enable us 
to decide if these pairs are true physical binaries. It is re- 
markable, though, that in both Lupus 3 and Chamaeleon I 
the binary candidates have a separation very close to the 
upper limit separation estimated for each cloud. 4 This 
might indicate that there is no real physical connection 
between the two members of these visual pairs. 

5. Ha emission 

As explained in Sect. |31 the value of the Ha-R colour in- 
dex is related to the equivalent width of the Ha line and 
hence, to the intensity of the Ha emission. The distribu- 
tion of the Ha-R colour with the spectral type is shown 
in Fig. The data seem to indicate an increase of the 
Ha emission with later spectral type, but the lack of ob- 
jects of lateM spectral type does not allow a definitive 
conclusion. Even if this were the case, we caution that 
the apparent higher emission might just be the effect of 
a lower continuum. Moreover, the two possible L-type ob- 
jects have very different values of the Ha-R colour index, 
which could be an indication of contamination with earlier 

4 Note that in Chamaeleon I this cutoff was larger, 12", due 
to its lower object density. 
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Fig. 7. Measured Hq-R colour index versus spectral type for 
our new Lupus 3 members. Spectral types 4 to 9 correspond 
to M4 to M9, while 10 to 12 correspond to L0 to L2. 
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M Subspectral Type 



Fig. 8. Upper panel: Measured Ha pseudoflux versus spec- 
tral type for our new Lupus 3 members. Lower panel Ratio 
LHa/Lboi versus spectral type. In both panels, spectral types 
4 to 9 correspond to M4 to M9, while 10 to 12 correspond to 
L0 to L2. 



type, highly extincted objects. Even with these uncertain- 
ties, it is remarkable that no decrease of the emission (and 
thus probably also the accretion rate, as discussed below) 
is seen across the substellar boundary. 

This can be better analysed by directly considering 
the Ha flux and/or luminosity of our objects. Following 
|Paper l| we computed a pseudoflux Fh q as: 



Fhq = Fhoj/Fq — 10 



_ i n-mn a /2.5 



(5) 



where mua denotes the Ha (instrumental) magnitude, and 
Fo is a hypothetical Ha absolute flux. We remark that this 
is a flux derived from the magnitude measured in our Ha 
band, and not the flux in the Ha spectral emission line, 
which can only be determined with spectroscopy. From 
this quantity Fh q , a Ha "pseudoluminosity" Lna was com- 
puted using the distance estimate of 140 pc. 

The upper panel in Fig. [S] shows the Ha pseudoflux 
as a function of the spectral type for our Lupus objects. 
The data apparently indicate a decrease of log Fh q with 
later spectral type. This is probably just the expression of 
the progressive faintncss of our objects. Indeed, if we plot 
LHa/Lboi instead, which should not be dependent on the 
distance or radius of the objects, no significant decrease 
of the emission with the spectral type is seen any longer 
(see lower panel in Fig. (SJ. An analogous behaviour was 
observed also in the Chamaeleon I cloud (cf. Fig. 14 in 
Paper 1). To obtain the bolometric luminosities, we have 
followed the method outlined by Comeron et al. (|2)00), 
using a distance value of 140 pc and R-I intrinsic colours 
from Kenyon & Hartmann l|1995f) and Zapatero Osorio et 

ai. C223. 

Since the Lupus members have Ha emission proper- 
ties similar to the Chamaeleon I objects, this emission 
must have a similar origin in both clouds. In the case of 
Chamaeleon I, the strongest emitters were shown to have 
a reported mid-infrared excess; hence, their Ha emission 
is most likely due to accretion from a circum(sub)stcllar 
disk. It seems thus plausible that also the Lupus objects 
are accreting from a disk. It is remarkable in this con- 
text that only one of our objects detected by 2MASS 
may have a near-infrared excess, a common indication of 
the presence of a disk around a pre-main sequence star. 
Nonetheless, several recent works (e.g. Comeron et al. 
l2rW)l Natta & Testi lWTTl Natta et al. l2nn2l|PapeTfl show 
that clear near-infrared excess is not always detected in 
objects with such cool photospheres. Hence, observations 
in the mid- and far-infrared might reveal disks around the 
Lupus objects. 

If confirmed as an intrinsic property of the objects 
(and not caused by a wrong distance or age estimation), a 
further indication of ongoing accretion would come from 
the apparent underluminosity of our new Lupus members 
(Sect. 1X2)1 . Undcrluminous objects have been reported by 
other authors (e.g. Fernandez & Comeron 120011 Luhman 
et al. I2TM1 Comeron et al. |211U5|) . They typically show 
indications of accretion (and sometimes also mass loss) 
in their spectra, particularly a prominent Ha emission 
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line. Comeron et al. H'2003|) discuss several possible origins 
for this apparent underluminosity, such as light scattering 
by a circum(sub)stellar disk seen edge-on or extinction 
through a dusty envelope. They finally suggest that the 
evolution of these objects may be critically affected by the 
accretion process, thus making them appear fainter than 
the predictions of the theoretical models, which are based 
on too simplified assumptions. Clearly, spectroscopy of our 
Lupus objects will be very useful to clarify their nature. 

Outflows or chromospheric magnetic activity may also 
contribute to the observed amount of Ha emission from 
our objects. It is not possible, however, to disentangle the 
mechanisms possibly responsible for this emission with the 
available data. 

6. Conclusions 

We have performed a survey for very low-mass objects in 
the Lupus 3 dark cloud. In an area of about 1.6 0° we 
found 19 new stellar and 3 brown dwarf candidate mem- 
bers. 

Like in Chamaeleon I, the Lupus low-mass objects are 
mostly located near the cloud core, with no evident signs 
of luminosity segregation. There is no obvious decrease of 
the Ha emission for the objects with the latest spectral 
types. The emission properties of the new Lupus members 
resemble those of the Chamaeleon I members. This indi- 
cates a similar origin, most probably accretion. However, 
in our analysis of the 2MASS data we find only one ob- 
ject which could have a near-infrared excess. Nonetheless, 
given that circum(sub)stellar disks are not always revealed 
by near-infrared photometry, we expect that many of our 
objects with Ha emission (if not all of them) may turn 
out to have disks in future mid- and far-infrared surveys. 

The most important caveats in our study are, on one 
side, the uncertainties in the distance and extinction of the 
Lupus 3 members, and, on the other, the uncertainties in 
our photometric spectral type classification, especially for 
the faintest objects in our sample, that might be affected 
by extinction. We particularly remark that in general our 
objects appear too faint for their estimated spectral type 
according to the theoretical models. While the most prob- 
able reason for this discrepancy is a wrong distance es- 
timate, we cannot rule out that the underluminosity is 
intrinsic to the objects and caused by the presence of cir- 
cumstellar disks seen edge-on. It is thus mandatory to 
perform spectroscopy to better constrain the true nature 
of our candidates and to further study their properties. 
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